Topological solitons in highly anisotropic two dimensional ferromagnets 
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We study the solitons, stabilized by spin precession in a classical two-dimensional lattice model 
of Heisenberg ferromagnets with non-small easy-axis anisotropy. The properties of such solitons 
are treated both analytically using the continuous model including higher then second powers of 
magnetization gradients, and numerically for a discrete set of the spins on a square lattice. The 
dependence of the soliton energy E on the number of spin deviations (bound magnons) A'' is calcu- 
lated. We have shown that the topological solitons are stable if the number A'^ exceeds some critical 
value A'cr- For A'^ < A'cr and the intermediate values of anisotropy constant Kcs < 0.35 J (J is an 
exchange constant), the soliton properties are similar to those for continuous model; for example, 
soliton energy is increasing and the precession frequency a; (TV) is decreasing monotonously with 
growth. For high enough anisotropy /Teff > 0.6J we found some fundamentally new soliton features 
absent for continuous models incorporating even the higher powers of magnetization gradients. For 
high anisotropy, the dependence of soliton energy E{N) on the number of bound magnons become 
non-monotonic, with the minima at some "magic" numbers of bound magnons. Soliton frequency 
iL>{N) have quite irregular behavior with step-like jumps and negative values of uj for some regions 
of A^. Near these regions, stable static soliton states, stabilized by the lattice effects, exist. 
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I. INTRODUCTION 



An analysis of two-dimensional (2D) magnetic soli- 
tons has been an a ctive area of research for more then 
30 years, see Refs. 11120 31 14151 for review. Such solitons 
are known to play an important role in the physics of 
2D magnetic systems. In the easy-plane magnets with 
continuously degenerate ground state there appear mag- 
netic vortices, responsible for the Berezinskii-Kosterlitz- 
Thouless (BKT) phase transitionSii. The presence of 
vortices leads to the emergence of a central peak in 
dynamical response functions of a magnet)^ which can 
be observed experimentallj*^. There are no vortices in 
the easy-axis magnets with discretely degenerate ground 
state, but various types of localized topological solitons 
appear there. Belavin and Polyakov were the first to con- 
struct the exact analytical solutions for 2D topological 
solitons in a continuous model of the isotropic magnetic. 
The energy of such Belavin - Polyakov (BP) solitons Ebp 
in a magnet with exchange constant J is finite and de- 
scribed by the universal relation 



^BP = 4:TTJS^, 



(1) 



S is the atomic spin. The structure of such solitons is de- 
scribed by a topologically nontrivial distribution of the 
magnetization field rh(x, y)^ which is determined by the 



TT2 topological invariant, see below Eq. H12|) and RefsHTI 
for more details. Belavin and Polyakov have also proved 
that such solitons are responsible for the destruction 
of the long-range magnetic order in purely continuous 
isotropic models at any finite temperature^ The studies 
of topological solitons has become interesting now due to 
their possible application in high-energy physics^ and 
the quantum Hall effecti^. Also, such soliton solutions 
after the change y — > ct (r is an imaginary time, c is the 
speed of magnons) determine so-called instantons, de- 
scribing non-small quantum fluctuations in ID isotropic 
antiferromagnetsi^. 

Note that the properties of the solitons in this isotropic 
model is rather academic problem since all real magnets 
have a discrete lattice structure and non-zero anisotropy. 
The role of uniaxial anisotropy has been investigated in 
a number of articles, see. Refs. Q|3 for review. Topo- 
logical solitons of the same 7r2 topological structure are 
also inherent for standard continuum (with accounting 
of terms quadratic on the magnetization field gradients, 
like the term W2 in the first line of Eq. ^ below) mod- 
els of anisotropic magnets. The basic problem of soli- 
ton physics in 2D magnets is related to soliton stability. 
According to the famous Hobard-Derrick theoremfi^*i& 
the stable static non-one-dimensional soliton with finite 
energy and finite radius does not exist in the standard 
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nonlinear field models; the soliton is unstable against 
collapse. This is, in particular, true for the uniax- 
ial 2D fcrromagnct with the anisotropy energy density 

The possibility to construct non-one-dimensional soli- 
tons stable against collapse is due to the presence of ad- 
ditional integrals of motion. For example, such solitons 
can be realized in the uniaxial ferromagnet due to the 
conservation of z-projection of the total spini*^. This 
leads to the appearance of so-called precessional solitons 
characterizing by time-independent projection of magne- 
tization onto the easy axis (z-axis hereafter), with the 
precession of the magnetization vector rh at constant fre- 
quency around the z axis. The analogies of such proces- 
sional solitons are known to occur in different models of 
field theory and condensed matter physics, see Ref. 
for review. Here the stability of 2D solitons is not re- 
lated directly to their topological properties; the non- 
topological dynamical solitons may also exist in magnets. 
Such solitons are characterized by the spatially localized 
function m(x, y), which does not have nontrivial topolog- 
ical properties. All stable precessional solitons, topolog- 
ical and non-topological, realize the minimum of energy 
for a given number of spin Sz deviations. In the semiclas- 
sical approximation, this value can take integer values N 
only, and can be interpreted as a number of magnons 
excited in the magnet. Thus, we naturally arrive at the 
concept of a soliton as a bound state of a large number 
of magnonsi. 

Two - parameter (parameters are precession frequency 
and velocity of translational motion of a soliton) small- 
amplitude non-topological magnetic solitons moving with 
arbitrary velocity in a two-dimensional easy - axis ferro- 
magnet have been constructed in RefllSi Their mini- 
mal energy -Ent = 11.7JS'^ depends only on combina- 
tion JS^ similar to Eq. which is a bit smaller then 
the energy of above discussed Belavin-Polyakov soliton, 
i?NT = 0.93i?BP- For such solitons, the relation between 
their energy (for given N) and momentum can be thought 
of as their dispersion law. Near the minimal energy 
i?NT, this dispersion law has a form E{P, N) ~ £{P) ■ N, 
where e{P) is a dispersion law for linear magnons. For 
P ~ 0, this gives the critical number of bound magnons 
Nc = £'NT/e(0). Thus, these solitons are nothing but 
weakly coupled magnon clouds, see Fig. [21 below. The 
expression for the soliton dispersion law is used to calcu- 
late the soliton density and the soliton contributions to 
thermodynamic quantities (response functions) like spe- 
cific heat. The signature of soliton contribution to the 
response functions of a magnet is an Arrhenius temper- 
ature dependence like exp^^Eo/T) with the character- 
istic value Eq as a soliton energy. Such behavior with 
i?NT < E(] < Ebp has been observed experimentally in 
Refs. I19I20I see Ref. for review. Comparison of contri- 
butions from solitons and free magnons shows that there 
is a wide temperature range where the solitons give more 
important contribution to thermodynamic functions such 
as heat capacity or density of spin deviations. 



We note here, that the structure of these non- 
topological solitons for E > E'nt is essentially differ- 
ent from that of topological solitons in uniaxial mag- 
nets. Namely, as ^ i^BP, the radius of topologi- 
cal soliton in continuous model of anisotropic magnets 
diminishes, making them "more localized", contrary to 
non-topological solitons, which become delocalized as 
E -Ent- Hence, although the energy of topological 
solitons is a little larger then that of non-topological 
solitons (0.89 and 1 in the units of Ebp, see above), it 
is possible that only BP-type topological solitons would 
contribute to response functions measured by neutron 
scattering in the region of non-small momentum trans- 
fer. On the other hand, recent Monte-Carlo simulations 
for 2D discrete models of easy-axis magnets did not show 
any signatures of small-radius topological solitonsSi. 

For real magnets, which are discrete spin systems on 
a lattice, there is an additional problem of application 
of the topological arguments, which, strictly speaking, 
can be applied only to continuous functions m(r, t). It is 
widely accepted that the continuous description is valid 
for discrete systems if the characteristic scale Iq of mag- 
netization rn(r,t) variation, |Vm| ~ m//o, is much larger 
than lattice constant a. The analysis of the magnetic vor- 
tices have shown that tti— topological charge of a vortex 
is determined by the behavior at infinities only so that 
the topological structure of such a vortex survives even 
in "very discrete" models with Iq < a. As for our case 
of 7r2 topological charge, the situation is not so simple 
and obvious. From one side, the continuous approach 
describes quantitatively the magnetization distribution 
in the vortex core already at Iq « I.SqSS. On the other 
side, so-called cone state vortices, with different energies 
for two possible spin directions in the vortex core are 
much more sensitive to the anisotropy, in fact, to the 
parameter o/Iq. Even for Iq > 10a their ^2— topological 
charge (polarization p = mz(0) = ±1), characterizing the 
core structure of vortices, for heavy vortices with higher 
energy can change so that they convert into more prefer- 
able light vortices with opposite polarization that 
never happened for continuum modc^i. Thus, the role 
of the discreteness effects is quite ambiguous. 

The above situation resembles the one-dimensional 
case - there are also kinks and breathers (non-topological 
solitons) with smaller energy. But it is well-known, that 
in the ID case for low anisotropy only kinks contribute 
to the response functions since small energy breathers 
transit continuously into weakly coupled magnon con- 
glomerates. We note that for such systems in one 
space dimension, the difference between topological and 
non-topological solitons for high anisotropy is not that 
large. For instance, the spin complexes with several 
A'^ ^ 10 magnons have been observed in a chain ma- 
terial CoCl2- 2H2O with high Ising - type anisotropjiS^. 
These complexes can be interpreted as non-topological 
one-dimensional solitons. 

The present work is devoted to the analysis of 2D 
solitons, both topological and non-topological, in the 
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strongly anisotropic magnets accounting for discreteness 
effects. In other words, here we investigate the influence 
of finiteness of a/Zo on the soliton structure. For inter- 
mediate values of anisotropy we found the presence of 
the critical number of bound magnons A^cr^ the topo- 
logical soliton is stable at > A'cr only. For very large 
anisotropy we found the specific effects of non-monotonic 
dependence of soliton properties on the number of bound 
magnons, caused by discreteness leading to the presence 
of "magic" magnon numbers. 



II. THE DISCRETE MODEL AND ITS 
CONTINUOUS DESCRIPTION 

We consider the model of a classical 2D ferromag- 
net with uniaxial anisotropy, described by the following 
Hamiltonian 



(2) 



Here S = {S'^, S^, S^) is a classical spin vector with fixed 
length S on the site n of a 2D square lattice. The first 
summation runs over all nearest-neighbors a, J > 
is the exchange integral, and the constant k describes 
the anisotropy of spin interaction. In subsequent dis- 
cussion, we will refer to this type of anisotropy as ex- 
change anisotropy (ExA). Additionally, we took into ac- 
count single-ion anisotropy (SIA) with constant K. We 
consider the z-axis as the easy magnetization direction 
so that K > or K > 0. 

In quantum case, the Hamiltonian (jSJ commutes with 
z— projection of total spin. It is more convenient to use 
semiclassical terminology, and to present it as a number 
of bound magnons in a soliton N, defined by the equation 



iV - ^ (5 - 55) 



(3) 



The spin dynamics is described by Landau-Lifshitz 
equations 



dt 



9Jf 

dSrl 



(4) 



the magnetization varies slowly in a space. In this case 
we can introduce the smooth function S{x,y,t) instead 
of variable Sn (t) and use a continuous approximation for 
the Hamiltonian It is based on the expansion of a 
classical magnetic energy E in power series of magneti- 
zation S gradients, 

E =.W2 + W4 + (5) 
where W2 contains zeroth and second order contributions 
to magnetic energy and W4 contains the fourth powers. 
These are given by 



W2 



Wa 




(6a) 



(6b) 



where V is a 2D gradient of the function S{r,t). Here, 
we used integrations by parts with respect to the fact 
that our soliton texture is spatially localized. We omit 
unimportant constants and limit ourselves to the terms 
of fourth order only as they are playing a decisive role 
in stabilization of solitons, see for detailaSiS^. Also, we 
introduce the effective anisotropy constant 



KeS = K + 2k. 



(7) 



We note here, that single-ion anisotropy enters only W2 , 
but not Wi and higher terms, while exchange anisotropy 
enters every term of the above expansion over powers of 
magnetization gradients. Usually, this difference is not 
important for small anisotropy, k ^ J, but, as we will 
see below, this fact gives qualitatively different behavior 
near the soliton stability threshold. 

Introducing the angular variables for normalized mag- 
netization 



TO = — = (sin 9 cos (f>; sin sin (/>; cos ( 



(8) 



In the case of weak anisotropy, K, k <ti J, the character- 
istic size of excitations Iq ^ a, see Eq. H14|) below, so that 



we obtain the following form of the classical magnetic 
energy 
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E\6, 



{yef (1 + K sin^ e) + iytpf sm^ ■ 



+ sin^ e (V(y9)^ [(V<^)^ + 2 (V6i)^j + 2 sin 6* cos 6{\I^6) \k (Vef - (Vip) 



(9) 



r 



In this long expression we omitted the terms with scalar 
product of gradients like (V6', V(/)), because they do not 
contribute to the trial function we will use for analysis, 
see below. 

In terms of fields 6 and cj), the continuous analog of 
Landau-Lifshitz equations (0J) read 



smi 



lyV^ a" SE „iyu 

m^hsJe' (10) 



^09 
dt 



a- SE 

Ks' 



These equations can be derived from the Lagrangian 

£[0,0]=-^ Jd'x{l-cos9)^-E[e,^]. (11) 

The simplest nonlinear excitation of the model pi|) 
with E = W2 is the 2D localized soliton, characterized by 
the homogeneous distribution of magnetization far from 
its core. Topological properties of the soliton are de- 
termined by the mapping of physical XY-plane to the 
§^~sphcre given by the equation = 1 of the order 
parameter space. This mapping is described by the ho- 
motopy group 7r2(S^) = Z, see^ii which is characterized 
by the topological invariant (Pontryagin index) 

Q = ^ y d^.T eai3 [m ■ (Vam x V/jto)] = 

= ^ J smO d6 d<f), (12) 

taking integer values, Q G Z. Here Sap is Levi-Civita 
tensor. 

To visualize the structure of a topological soliton, 
we consider the case of a purely isotropic magnet with 
Kcs = and W4 = 0. For this isotropic continuous 
model the soliton solution is aforementioned BP soliton 
of the formifi 



tan ■ 



'Po + Qx, 



(13) 



where r and x ^re polar coordinates in the XY-plane, ipQ 
is an arbitrary constant. The energy of this soliton 
does not depend of its radius R, which is arbitrary pa- 
rameter for the isotropic magnet, see Ref. ITol However, 



even small anisotropy breaks the scale invariance of the 
above model since now the scale 



2K 



cff 



(14) 



enters the problem. 

In the latter case, the soliton energy with respect to 
W2 term only has the form E = Ebp + const • (R/Iq)'^ 
and hence has a minimum at i? = only, which signifies 
the instability of the static soliton against collapse in this 
model. Obviously, the scale invariance is also broken in 
the initial discrete model. The "trace of discreteness" in 
our continuous model Hll|l is the presence of a contribu- 
tion W4. This term gives the contribution to the soliton 
energy proportional to JS^a^ /R^. For some rather exotic 
models containing higher powers of a magnetization field 

gradients like (^nij (with positive sign) (22i224^2i^ such 

term might be able to stabilize even static soliton against 
collapse, seoS^ for details. However, for Hcisenbcrg mag- 
nets this kind of stabilization is very problematic. For 
example, discrete magnetic models with the Heisenberg 
interaction of nearest neighbors only (i.e. those without 
biquadratic exchange and/or next-nearest-neighbors in- 
teraction) have negative W4 < 0, and the higher powers 
of magnetization gradients do not stabilize a static soli- 
ton in this case. Moreover, as we will show below, the 
presence of discreteness ruins the stability of the preces- 
sional soliton with N < N^- even in the case, when it is 
stable in the simplest model with W = W2 only. 

We now discuss the stability of precessional solitons. 
For a purely uniaxial ferromagnet the energy functional 
E[9,(j)] does not depend exphcitly on the variable so 
that there exists an additional integral of motion^i which 
is the continuous analog of Eq. ^ 



N 



S 



d'^x{l~ cosd) . 



(15) 



The conservation law (|15|l can provide a conditional (for 
constant N) minimum of the energy functional E, stabi- 
lizing the possible soliton solution. Namely, we may look 
for an cxtremum of the expression 



(16) 



where lo is an internal soliton precession frequency, that 
can be regarded as a Lagrange multiplier. Note that this 
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functional is nothing but the Lagrangian Hll() calculated 
with respect to specific time dependence 

(t) = ujt + Qx + fo, 
which holds instead of (|13|l in this case. This condition 



hiu = 



dE 
dN' 



(17) 



which makes clear the microscopic origin of the preces- 
sional frequency uj. Namely, an addition of one extra 
spin deviation (bound magnon) to the soliton changes 
its energy by fku. Thus, the dependence u^N) is ex- 
tremely important for the problem of a soliton stability. 
For the general continuous model of a ferromagnet, even 
containing the terms like W4, the sufficient and necessary 
condition of soliton stability reads dui/dN < 0^ but for 
the discrete model the validity of this condition is not 
clear yet. The point is that the known analytical meth- 
ods of a soliton stability analysis rely essentially on the 
presence of a zcrotli (translational) mode, which is obvi- 
ously present for any continuous model, but is absent for 
discrete models, where lattice pinning effects are present. 
To solve the problem of soliton stability we will investi- 
gate explicitly the character of conditional cxtrcmum of 
the energy with N fixed. 



III. THE METHODS OF A SOLITON 
STRUCTURE INVESTIGATIONS 

To get the explicit soliton solution and investigate its 
stability wc have to solve the Landau-Lifshitz equations 
(|10|l with respect to the energy l^. For the simplest 
model accounting for W2 only, an exact ansatz 



(r), (l3 = Qx + ujt , 



(18) 



can be used, leading to the ordinary differential equation 
for the function 0{r). This equation can be easily solved 
numerically by shooting procedure, using the value of 
dO/dr at r = as a shooting parameter, sec Ref. Q. The 
shape of the soliton essentially depends on the number TV 
of bound magnons. In the case of the soliton with large 
N ^ N2 = 2tt S {Iq / a)'^) , the approximate "domain wall" 
solution works pretty well. This solution has the shape 
of a circular domain wall of radius R 



cos 6*0 (r) — tanh 



r- R 



(19) 



Using this simple structure one can obtain the num- 
ber of bound magnons, which is proportional to the 
size of the soliton, « 2TTS(R/af and energy, E = 
47TS^^/2JKcsR/a. Note that such a solution is the same 
for topological and non-topological solitons except for the 
behavior near the soliton core. This means, that the 
characteristics of such solitons are pretty similar. In the 



case of the small radius soliton {R <C Iq), the following 
asymptotically exact solution works well^i 



tan = — A 1 — 

2 ro V ''0 



^0 



^1 - LO/UJQ 



(20) 



where Ki(x) is the McDonald function, and loq is a gap 
frequency for linear magnons. It provides correct behav- 
ior at r < i? ^ Zo> where it converts to the Belavin- 
Polyakov solution (|13|l . For large distances {r ^ R), this 
expression gives an exponential decay (instead of power 
decay for a Belavin-Polyakov soliton) with characteris- 
tic scale tq. For solution (|20|l u; ^ ujq as N ^ so that 
small radius solitons in anisotropic magnets have two dif- 
ferent scales, the core size R Iq an d the scale of the 
exponential "tail" ro = Iq / ^/T^~LuJLJ^ > l^^i^. 

For even minimal accounting for the discreteness on the 
basis of a generalized model with fourth spatial deriva- 
tives, the problem becomes much more complicated. The 
complexity of the problem is not only due to the fact 
that for the energy © it is necessary to solve the fourth 
order differential equation, and use a much more com- 
plicated three-parameter shooting method. But the ba- 
sic complication here is the fact that in general fourth- 
derivative terms contain anisotropic contributions like 
{d'^O/dx^Y + {d^d/dy'^Y, which cannot be reduced to 
the powers of radially-symmetric Laplace operator (these 
terms are omitted in Eq. for simplicity). In this case, 
the radially-symmetric ansatz with 6 — 0{t), (j) = 4>{x) is 
not valid, and we have to solve a set of partial nonlinear 
differential equation for the functions 9 = 9{r,x), 4' = 
4>{r, x)- To the best of our knowledge, an exact method 
for construction of soliton (separatrix) solutions for such 
type of equations do not exist so that some other approx- 
imate methods should be used for this problem. 



Variational approach for general continuum 
model 



One of the approaches, which we use for the approx- 
imate analysis of the solitons in the model is the 
direct variational method. For the minimization of the 
energy E = 11^2 -I- we use a trial function. 



tan - = Ai?A:i (Ar) . 



(21) 



Here, we consider the case Q = 1, and the case of higher 
topological invariants is qualitatively similar. Note, that 
trial function 121|) is based on the asymptotically exact 
soliton solution (|20|l . The trial function H21|) gives cor- 
rect asymptotics both for x — > (corresponding to BP 
soliton) and for a; — *■ 00 (exponential decay with charac- 
teristic scale 1/A), see above. Our analysis shows that 
the same results can be obtained using a simplified trial 
function, which also captures the asymptotic behavior of 
the soliton. This function has the form 



tan ■ 



R 



exp(— Ar) 



(22) 
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In the spirit of the minimization method above discussed 
we consider the parameter A as variational, keeping the 
parameter R constant as it is related to A^, N cx R^, 
see, e.g.i In other words, we minimize the energy (PJ 
with the trial function (|21|l or (|22|l over A for constant 
R. This approach has the advantage that simultaneously 
with equation solving, it permits investigation of the sta- 
bility of obtained solution on the base of simple and obvi- 
ous criterion. Namely, a soliton is stable if it corresponds 
to the conditional minimum of the energy at fixed A'^, and 
it is unstable otherwise. 



Numerical analysis of the lattice model 
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Since the continuous description fails for the case of a 
high anisotropy, one needs to elaborate the discrete en- 
ergy (O on a lattice. As we are not able to solve the 
problem analytically, we discuss the numerical approach 
to study the soliton-like spin configurations on a discrete 
square lattice. Obviously, the direct molecular dynamics 
simulations on a lattice is a powerful tool for the soliton 
investigation for any values of anisotropy. Direct spin 
dynamics simulations of 2D solitons have been recently 
performed in Refs. Is^ls^ . After the relaxation scheme, 
where an initial trial solution was fitted to the lattice, 
the spin-dynamical simulations of the discrete Landau- 
Lifshitz equations were performed. This approach, how- 
ever, is quite computer intensive, requiring powerful com- 
puters. To minimize the calculation time, the parallel 
algorithms have been used24. 

For this reason, we obtained desired spin configuration 
by direct minimization of the energy E, keeping A'^ con- 
stant. For our modeling we choose the square lattices 
with "circular" boundary conditions. We fix the values 
of spins on the boundary to the ground state {9 = 0), 
those values have been kept intact during minimization. 
The method of minimization is the simplex type method 
with non-linear constraints. This method is based on 
the steepest descent routine applied to the functions of 
a large number of variables. The above method is able 
to find the conditional minimum of a given function with 
several (usually small number) constraints, consisting of 
relations between variables. In our problem, such vari- 
ables are the directions of each spin, parameterized by the 
angular variables 0„ and (/3„, and the conditional mini- 
mum of the energy, E has been obtained for a fixed value 
of z— projection of the total spin. As this additional con- 
straint slows down the calculations substantially, we use 
other method, valid for initial configurations where A" 
differs from the necessary quantity by one spin. The idea 
of the method is as follows: the angle 9 of an arbitrarily 
chosen "damper" spin was excluded from the minimiza- 
tion procedure, and its value has been kept constant, to 
achieve the necessary A^ value throughout whole mini- 
mization with respect to all other variables. 

The method in fact is dealing formally with a static 
problem, but it gives the possibility to find the preces- 



FIG. 1: Dependence E{N) and u){N) for a soliton with KcS = 
0.3, obtained from numerical simulations. 



sional frequency directly. To find uj, we used the discrete 
Landau-Lifshitz equation rewritten for the angular 
variable in the way used in equation (|10|l . Using 
relation dcjjn/dt = to, we obtain 

- ^ 1 d:K 
~ sin6'„ ' hSdK ' 

which docs not depend on index n throughout all the 
system. Here J{ is the discrete Hamiltonian |j2J). 

To find the above local minimum, we start from the 
BP initial configuration. The size of lattice clusters var- 
ied from 20x20 (for large anisotropies, where discreteness 
of a lattice is revealed most vividly) to 32 x 32 for small 
anisotropy when system is well described by a continuous 
model. The criterion of the presence of a truly local soli- 
ton configuration was its independence of the system size. 
Another important criterion was the constancy of the fre- 
quency w, calculated from the equation H23() throughout 
the soliton configuration. Sometimes we found the min- 
imum over the variables 9n only, considering 0„ to obey 
Eq. (|18|l and choosing the reference frame origin in the 
symmetric points between the lattice sites. Our analysis 
has shown that for moderate anisotropies K^s ^ 0.5J 
the above partial minimization gives the same energy 
and frequency as well as the instability point position, as 
complete minimization over 0„ and . The use of par- 
tial minimization over 9n permits not only to accelerate 
the numerical calculations, but turns out to be useful for 
construction of " quasitopological" textures for extremely 
high anisotropies, see below. Section^] In FigQ] the de- 
pendence of a soliton energy and precession frequency, 
is shown as a function of the bound magnon number for 
A'cft — 0.3 for single-ion and exchange anisotropies. The 
frequency has been calculated directly from the equation 
(|23|l and also found by differentiation of the energy with 
respect to A^, see Eq. H17(l . Here, far from the instabil- 
ity point A^cri the behavior of uj{N) is almost similar, 
while near this point it is different. For both types of 
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anisotropy (single-ion and exchange) the "cusps" in the 
dependencies uj{N) a.s N ^ N^r are well seen. 

For N < Nci the minimal configuration cannot be 
found. In fact, any attempt to find the topological soli- 
ton with N < Ncr leads to appearance of non-topological 
solitons with dip/dx — 0. This occurs even for the above 
"partial" (i.e. over 0„ only) minimization. In this case, 
some spins rotate by non-small angles to organize the 
structure with zero Q. This demonstrates the instability 
of topological solitons for TV" < A^^r quite vividly. These 
features, as well as the values of A^cr or E^, can be de- 
scribed also on the basis of a variational approach with 
simple trial functions, see the next section. 

IV. SOLITONS FOR MODERATE VALUES OF 
ANISOTROPY. 

To perform specific calculations for the generalized 
continuous model H6a|) . it is convenient to introduce the 
following dimensionless variables 

X = Ar, A = aA, z = Ai?, (24) 



as before, a is a lattice constant. Further we may ex- 
press both trial functions (|21|l and (|22|l in the following 
universal form 

tan^ = zf{x), h{x) = K^{x), h{x) = ^^Ziz^. (25) 
2 X 

Then using the equations (|15|l and ©, we can calculate 
the number of bound magnons in the soliton 

- = iir-^Piz), yj{z) = / xdx J ^ > , (26) 

and the soliton energy 

_|_.K£,.m+,,„__LaSW^ (27) 

Here we introduced the following notations 



' . 1 f°° 

sin^ 6xdx, 12(2) — — / xdx 
2 Jo 



(l + Ksin^i 



sin^6l 



74(2) = J xdxl {A^ef (1 + Ksin^ 6) + 9"^ (l + kcos^ 9) 



sm 



+ 29'- 

x'^ 



+AJsm29 k9' 



1 



de 



d^9 1 d9 



dx ' ^ dx^ ' X dx 



(28) 



Thus, we express the energy and the number of 
magnons via two parameters, A and z. It turns out, that 
initial dimensional variables A and R enter the problem 
only in the form of their product z. The dependence of 
A^ and E on z enters the problem via a few complicated 
functions ■0, 70i 72 and 74, which can be written only im- 
plicitly in the form of integrals. However, in terms of 
these functions, the dependence on A (|27() turns out to 
be quite simple. This permits reformulation of the initial 
variational problem in terms of variables z and A'^ only. 
Namely, we express 

A^=4.^^(.) (29) 

and substitute this expression in the dimensionless en- 
ergy l|77|l . This gives us the expression for the energy of a 
soliton with given A'^, as a function of variational param- 
eter z. Then we can find a minimum of E with respect to 
z, keeping A^ constant. The result of such minimization 
in the form of the dependence i?inin(A^) is shown in Fig[21 



for a magnet with purely exchange anisotropy (i.e. for 
that with K = 0). We show this dependence for k = 0.2. 
It is seen that there is a good correspondence between the 
dependencies E{N) found by variational and numerical 
(full symbol) minimizations in the region of parameters 
where the soliton is stable. This justifies the applicabil- 
ity of the variational approach with the trial functions of 
the form H25|l to the problem under consideration. The 
approach based on simpler trial functions /2 reproduce 
well the particular feature found for numerical analysis 
of the discrete model; namely, the value of the threshold 
number A'cr- 

Also, for comparison, we show in Fig. [2 the result of 
variational minimization of W2 (i.e. energy, incorporat- 
ing only squares of magnetization gradients). It is seen, 
that in this case there is no A'cr, which coincides well with 
the previous investigations of solitons in the continuum 
modelsiSSiSi This means, that mapping of the initial dis- 
crete model on the simplest continuum model with only 
squares of magnetization gradients can be wrong for some 
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FIG. 2: Dependence Emin{N) for topological and non- 
topological solitons for the case of exchange anisotropy only, 
K = 0.2. Dashed lines, parallel to the axes, correspond to 
characteristic number of bound magnons N = Nc for non- 
topological soliton, see the text above. Curves, corresponding 
to contribution of W2 only are shown for comparison. 
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FIG. 3: Soliton phase diagram, (a) - exchange anisotropy, 
solitons exist in the region 1 and do not exist in the region 
2. For small k the value of A'cr ~ lAS/n. (b) - com- 
parison between exchange (curve 1) and single-ion (curve 
2) anisotropics. Limiting values of A'^cr for both types of 
anisotropics are shown by vertical dashed lines. 



values an N, and to get the correct description of solitons 
in a 2D magnet one must take into account at least fourth 
powers of gradients. On the other hand, for large enough 
N and even moderate value of anisotropy, the role of this 
higher derivative terms is less important; it is in agree- 
ment with the recent numerical simulations of soliton dy- 
namics for easy-axial discrete models of ferromagnets^. 

It is seen from Fig|21 that the topological solitons in 
such a model are not very "robust" - there is a quite 
large parameter region where those solitons arc unstable 
or do not exist at all. That is why for comparison we 
show the dependence E{N) for non-topological solitons 
(open symbols on the figure) which have been studied in 
detail analytically in Ref. The main feature of non- 
topological solitons is that while — > wq, the amplitude 
of such solitons diminishes (so that at lu > ujq such soli- 
ton decays to a number of noninteracting magnons) but 
both its energy E = i?NT ~ 11.7i?BP and bound magnons 
number N = Nc ~ E-^^/hwo arc left intact liS, This be- 
havior is opposite to that of a topological soliton, where 
as w — > cjQ the amplitude still has its maximal value with 
decreasing soliton radius so that soliton becomes "more 
localized" . This tendency is seen in Fig|21 where a non- 
topological soliton exists up to TV = N^t (corresponding 
to w = Wo, shown as vertical dotted line on the figure) 
and then decays smoothly into a noninteracting magnon 
cloud with the energy above on the corresponding curve 
(dotted line in the figure). This is because as uj ^ 
the non-topological soliton has the same energy as corre- 
sponding magnon cloud. But the soliton is " coherent" (in 
a sense that it is a bound state of many magnons), while 
magnon clouds are not coherent. Actually, the above re- 
mark refiects the important point in a physics of solitons 
under consideration, namely the difference in behavior of 
topological and non-topological solitons. 

As was mentioned above, there is a critical value of 
bound magnons number A^cr such that the topological 
soliton exists with N > iVcr only. This threshold de- 
pends on K. This dependence, which we will call a phase 
diagram, is depicted in FigOt. The solitons exist in the 
region 1 above corresponding curves and do not exist in 
the region 2 below them. For comparison, on the same 
figure, we plot the points iVcr, corresponding to differ- 
ent trial functions. Very good coincidence between these 
points shows that both trial functions are well-suited for 
variational treatment of the solitons. 

It is seen that the threshold number of bound magnons 
calculated using the simple variational approach grows 
infinitely both for k — > and for k — > Kum ~ 

0.312. 

The comparison with the numerical data shows that 
the divergence at the large values of anisotropy is sim- 
ply an artifact of continuous description, and the cor- 
responding "numerical" curve for discrete model de- 
creases monotonously with growing of anisotropy. For 
high anisotropies, the soliton structure becomes strongly 
anisotropic, and the description based on radially- 
symmetric trial functions fails. Wc will discuss this in 
the next section. 
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On the other hand, the divergence of Ncy at small k 
has a clear physical meaning and coincides well with the 
numerical data. The divergence of A^cr at k is re- 
lated to the fact that in BP soliton, which is the exact 
soliton solution for purely isotropic fcrromagnet, the in- 
tegral describing the value of bound magnons N diverges 
logarithmically as r — > oo due to slow decay of the func- 
tion cx 1/r. Using the variational approach with 
the asymptotics of the functions H28|l it is possible to de- 
rive the asymptotic formula for Ncr in this region, which 
reads 



iV„(K) « (30) 



In contrast, the critical value of energy Ea = E{Ncr) is 
finite at the instability point, 

E„{k) K 47rJS^{l + 2MyQl). (31) 

This limiting energy contains non-analytical dependence 
on the anisotropy constant k. It appears due to simulta- 
neous accounting of the anisotropy and fourth derivative 
terms in our generalized continuum model. For typical 
values k/J = 0.1 0.2, this energy is well above the 
Belavin-Polyakov limiting energy, see Figs. [5] The 
correction to Ebp becomes smaller then 1% at extremely 
low anisotropics like k < 10~^ J only. The situation here 
is very similar to the analysis of heavy (less energetically 
favorable) vortex decay in the cone state of an easy plane 
ferromagnet(2^ which situation occurs in a magnetic field, 
perpendicular to the easy planei2i24 There, in the sim- 
ple continuum model, the heavy vortices are stable in 
the entire region of cone state existence {0 < H < Ha, 
Ha is an anisotropy field)^, but already for very small 
anisotropics k ~ 10~^, when lo ~ 10a, this region has 
diminished substantially so that at k ~ 0.1 the heavy 
vortices have already become absent 

In Fig|2|3, the phase diagrams for exchange (curve 1) 
and single-ion anisotropics (curve 2) are shown. It is seen, 
that while at small k and K the corresponding curves 
lie close to each other, for larger anisotropics there is a 
drastic difference. While the (unphysical) limiting value 
of A^cr for high exchange anisotropy k equals 0.312, cor- 
responding to KcS = 0.624, the same value for single- 
ion anisotropy constant K corresponds to KcS ~ 0.5. 
This means that the continuous description for exchange 
anisotropy is valid for larger values of anisotropy con- 
stants, even in the region Kcs — 0.5, where it already 
fails for single-ion anisotropy. This is a consequence of 
the fact, that the single- ion anisotropy constant K enters 
the problem only in the spatially homogeneous term (via 
the combination Kes, see above), while the constant k 
enters all terms of expansion Wi . 



V. SOLITONS IN THE DISCRETE MODEL 
WITH HIGH ANISOTROPY. MAGIC NUMBERS 
OF BOUND MAGNONS 

For large anisotropics, the discreteness plays a deci- 
sive role so that exact analytical treatment of the prob- 
lem is impossible. However, an account for the fol- 
lowing two facts permits to accomplish a comprehen- 
sive approximate study. First, as was mentioned in the 
early publicationsr^i^^i^ for the large number of bound 
magnons > S(lo/a)'^ , any solitons in two- and three 
dimensional magnets, topological and non-topological, 
can be presented as a finite region of flipped spins, sep- 
arated by 180° domain wall (DW) from the rest of a 
magnet. 

Second, for such solitons (magnon droplets) the dif- 
ference between topological and non-topological textures 
becomes vague. As the anisotropy grows, Iq ^ a and the 
characteristic magnon number becomes comparable with 
2ttS. At the same time, the DW becomes thinner so that 
a.t K, K ^ J its width becomes comparable with lattice 
constant a. It is clear that already at high anisotropy 
the structure of bound states (solitons) with N/S > 10 
is above the described flipped area bordered by the DW 
(we recollect here that we consider a soliton as a bound 
state of many magnons) so that the soliton properties 
will be completely determined by those of the DW. It is 
also clear that for this case the difference between topo- 
logical and non-topological solitons will be negligible. We 
will see, that in this case the structure of bound state is 
strongly dependent on the character of anisotropy. That 
is why for a description of such bound states it is useful 
to study first the structure and properties of the DW's 
in highly anisotropic 2D magnets with different types of 
anisotropy. First of all, of importance is a notion of the 
DW pinning, i.e. the dependence of its energy both on 
the DW center position in a lattice (positional pinning) 
and on its orientation with respect to the lattice vectors 
(orientational pinning). 

The positional pinning of a DW can be discussed on 
the basis of a simple ID model of a spin chain. This 
case is obviously applicable to the 2D square lattice 
with nearest-neighbor interaction for the DW orientation 
along lattice diagonals (the directions of (1,1) type). For 
a chain, it is natural to associate the DW coordinate X 
with the total spin projection on the easy axis §tot°^ 
to define it as follows^. Let us choose some lattice site 
and define the DW located on this site to have the co- 
ordinate X = 0. Let us then determine the coordinate 
of any DW via its total spin z-projection Sf^j from the 
expression 

n—oo 

^ = ^ E [S^AX)~S^„{X^O)], (32) 

n— — oo 

where S^{X — 0) defines the spin distribution in a DW 
at a reference point X = 0. 

For the solitons description, the different properties 
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of DW placed in different positions play a crucial role. 
First, consider a DW in a middle between two arbitrary 
lattice sites, for which X/a is half-odd, X — a(2n+ l)/2. 
Only domain walls of such type in the highly anisotropic 
magnet can be purely coUinear, i.e. it may have Sz = ±5", 
respectively, on the left and right sides of the DW center. 
For the case of a spin chain with single-ion anisotropy 
this is achieved for K > 0.5p^ which was associated with 
the destruction of non-coUinear topological structure^. 
According to the definition l|32|) . it is a vital necessity 
to have a non-integer S^^j and non-coUinear component 
for the rest of the DW positions. For example, for X = 
(a/2)(2n) in the one of the sites Sz = 0, i.e. 9 = n/2. 

The question regarding the character of the pinning 
potential U(X) is also important. Gochev has shown, 
that there is no pinning, i.e. C/ = 0, in a spin chain with 
purely exchange anisotropy^ On the other hand, for the 
anisotropy of a pure single-ion type, the potential U{X) 
has minima in the points of type a{2n + l)/2, which fa- 
vors the appearance of coUinear DW'a^. For the small 
perturbation of a problem with exchange anisotropy by 
single-ion anisotropy with negative sign K < 0, it turns 
out that U (X) has the minimum at integer values of X/ a, 
X = 0, ±a, i.e. the creation of a non-coUinear struc- 
ture becomes favorable. The detailed investigation of the 
character of DW pinning for general models with com- 
bined anisotropy shows, hoverer, that this situation is 
mostly an exception and the intensity of such pinning is 
never large. The pinning potential for the parameters 
where it has minima at integer X/a is much smaller (at 
least one order of magnitude) than that for pure single- 
ion anisotropy. 

Thus, for the two above mentioned cases one can ex- 
pect substantially different spin distributions in a soliton 
boundary. In particular, DW pinning on lattice sites can 
yield the conservation, at least partial, of a soliton topo- 
logical structure even with strong discreteness effects. On 
the other hand, a DW pinned between lattice cites be- 
comes coUinear and hence its topological structure can 
be easily lost. 

It is clear, that for the question about the structure 
of a closed soliton boundary, the important point is its 
angular pinning, i.e. the dependence of the planar DW 
energy on its orientation in a crystal. Gochev suggested, 
that the optimal DW direction is parallel to the primi- 
tive vectors of lattice translation, ( 1,0) and (0,1);^ that 
differs from the conclusions of Ref. Our analysis has 
shown that almost all numerical data about the prop- 
erties of the solitons in highly anisotropic FM's can be 
described under supposition that a DW tends to be par- 
allel to (0,1) direction. As we shall show, only N/S is 
a main parameter determining the soliton structure at 
high magnetic anisotropy. It turns out that the variation 
of N/S by around few percent yields substantial varia- 
tions of DW structure, which leads to sharp dependence 
of energy E (see Fig^J and particularly the frequency 
u! on N, see Fig. [3 below. This effect is different for 
single-ion and exchange anisotropics, the most substan- 




N/S 



FIG. 4; Dependence E{N), obtained by numerical simula- 
tions for Kcff — 1. Curve 1- single-ion, K — 1, curve 2 - 
exchange, n — 0.5. 



tial manifestation is for the single-ion case. These de- 
pendencies have a non-monotonic component. Its analy- 
sis reveals certain specific numbers A^mag, which, analo- 
gously to nuclear physics, can be called the magic num- 
bers. To explain the origin of these magic numbers, we 
consider the case when a DW tends to occupy a posi- 
tion between atomic planes of (1,0) and/or (0,1) type so 
that both DW bend and non-coUinear structure forma- 
tion are unfavorable. Then, the optimal (from the point 
of view of DW energy) configuration, is that where the 
spins with Sz = —S occupy a rectangle Ix ■ ly, separated 
from the rest of a magnet by a coUinear DW. Clearly, 
the most favorable configuration is that with spin square 
Ix = lyj which yields Amagic = 2/^5, but configurations 
with Ix ~ ly; like those having N/S = 2 • + 1), for 
example. A/ 5 = 2 • (3 • 4) = 24 are also quite profitable, 
see Fig. 0] We will call such values of A^ half-magic. 
As we shall see, such a model describes the solitons in 
magnets with single-ion anisotropy for N/S ~ 15 50 
pretty well. If the energy of spatial and angular pinnings 
is not so important, we may expect more or less cylin- 
drical shape of a soliton core and smooth dependencies 
E{N) and uj{N). The numerical analysis has demon- 
strated, that both above tendencies reveal themselves in 
a FM with single-ion and exchange anisotropics respec- 
tively, see Fig. 21 for E{N) dependence. Namely, for SIA 
the dependence E{N) has a non-monotonous component, 
while for ExA this dependence seems to be more reg- 
ular. The minima in the non- monotonous dependence 
E{N) for SIA occur for N/S = 18 = 2 • 3^, 24 = 2(3 • 4) 
and 32 = 2 • 4^. Hence, the above magic and half-magic 
numbers, related to the DW pinning between two atomic 
planes, are clearly seen in the E{N) dependence. 

We now discuss non-topological spin configurations for 
different numbers of bound magnons for ferromagnets 
with single-ion and exchange anisotropy with large value 
of effective anisotropy constant ATgff = J, see FigsISjand 
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(a) ExA, 8 



(c) ExA, 12. 



(c) ExA, 14. 



(g) ExA, 18. 



(b) SIA, 8 
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(d) SIA, 12. 
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(f) SIA, 14 




(h) SIA, 18. 



FIG. 5: The structure of the sohton textures for the magnets 
with exchange (left column) and single-ion anisotropy (right 
column) with the same value of effective anisotropy and dif- 
ferent values of n = (shown as subfigure captions). The 
arrows present in-plane spin projections in 20 by 20 lattice. 
"Up" (0 < f < 10°) and "down" (170° < 6* < 180° ) spins 
are presented by dotted and crossed circles, respectively. The 
in-plane projections of the spins with "up" and "down" z — 
projections are depicted by arrows with open and full heads, 
respectively. 



Let US first consider the small values n = N/S < 20. 
To economize the notations, hereafter we will use n in- 
stead of N. For very small n < 10 the soliton textures 
have the same noncoUinear structure for both types of 

5(a)|5(b) 



The dif- 



anisotropies, SiA and ExA, see Figs, 
ference between spin textures for SIA and ExA solitons 
becomes visible for n > 10. In this case, the effects of 
lowering of a soliton symmetry with respect to expected 
lattice symmetry of fourth order C4, arc possible. For 
half-magic number n = 12, corresponding to collinear 
texture with flipped spins rectangle 2x3, the symmetry 
C4 is obvi ously absent both for SIA and for ExA, see Figs. 
5(c)|5(d)' Also, there is n o big difference between spin 
textures in Figs. |5(c)|5(d') However, for n far from magic 
numbers, the difference between SIA and ExA is much 
more pronounced. In particular, for the SIA case there is 

while for ExA, the C4 sym- 
The explanation is pretty 



not even a C2 axis (Fig. 
metry is restored (Fig. 
simple: the DW pinning is weaker for ExA then that for 
SIA so that in former case the symmetric closed DW is 
formed, while for SIA more favorable is the formation of 
a piece of "unfavorable" DW, occupying only the part of 
a soliton boundary, which makes possible optimization 
of DW structure for the rest of the boundary. As n in- 
creases further for SIA the purely collinear str uctures of 
the above discussed type, can appear, see Fig 5(h)|6(j) " 



As for the ExA case, even at sufficiently large n = 32, the 
soliton texture does not contain a purely collinear DW, 
see n = 18 in Fig 5(g) and large n = 32 in Fig 6(i) below. 



Thus, at small the certain tendency, which is con- 
firmed at large n (see Fig. ISJ, is clearly seen. Namely, the 
strong DW pinning for the SIA case yields almost always 
the nonsymmetric configurations, where n growth occurs 
due to increase or decrease of "DW pieces" on a soli- 
ton boundary. Figs 6(b)|6(f)] Th e exceptions are magic 



6(j) I or close t o them 
6(d) I, where 



(Fig. 



numbers rimag = 2Z^ (Figs, 
"half-magic" numbers rihm 
the collinear structure of the flipped spins square or rect- 
angle type (with symmetries, respectively, C4 and C2) is 
present in SIA case. For the ExA case the numbers nmag 
and rihrn are also revealed, but in a quite different man- 
ner. Namely, for n = nmae the soliton does not have a 



collinear structure but resembles a square, see Fig. 6(a) 
Here, however, there is a fundamental difference with 
the SIA case. The same structures with C4 symmetry 
are formed also for n non-magic, n ^ fT-mag, ^^hm- This 
becomes clear if we recollect that in this case the pin- 
ning is weak (see above) so that "DW piece" formation 
is absolutely unfavorable. This means that C4 symmetry 
occurs both for n magic and nonmagic, see Fig |5(g)|6(i) 



However, for n half-magic, the rectangular shape of a soli 
ton core occurs al so for ExA both for small (Fig 
and large n (Fig, 



6(c) I 



5(c) I 



Thus, the soliton structure for 
the ExA case also depends on n nonmonotonously, but 
here the "half magic" numbers are more important since 
close to these numbers the soliton symmetry first lowers 
from C4 to C2 and then it resto red b ack to C 4. At the 
magic number n = 32, see Figs. 6(i) and 6(j) as well as 



12 



(a) ExA, 21. 



(c) ExA, 24. 



(c) ExA, 28. 



(g) ExA, 31. 



(i) ExA, 32 
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(b) SIA, 21. 




(d) SIA, 24. 
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(f) SIA, 28. 




(h) SIA, 31. 




(j) SIA, 32. 



near this magic number, both SIA and ExA textures have 
quite symmetric structure with almost colhnear DW, see 
Figs. 6(g) and 6(h) for n = 31. 

This complex and irregular picture of soliton behavior 
"maps" onto the E{N) dependence only as a local energy 
lowering at n w ^T-mag- The irregular behavior of soliton 
characteristics is revealed much more vividly in the de- 
pendence uj{N), seen in Figs. |Zta),(b). For both types of 
anisotropy, the explicit traces of nonmonotonous behav- 
ior like jumps, regions with duj/dN > 0, and even those 
with w < (for SIA, see Fig[7J) occur. At this point it 
is useful to make two remarks. First, we recollect, that 
for discrete systems, contrary to the continuous ones, the 
condition du) /dN < is no more a soliton stability crite- 
rion. Second, the condition w < just means that in this 
region of parameter values, the soliton energy decreases 
as N increase, but says nothing about the soliton stabil- 
ity. It is seen from Fig Ufa), that the non-monotonous 
structure of the corresponding lj{N) curves manifests it- 
self most vividly atn > 10-^15. For SIA, the abrupt ver- 
tical up and downward jumps with growing of n on the 
curve uj{N) appear at certain values of n. The upward 
jumps occur near magic and half-magic numbers intro- 
duced above. After these upward jumps, the frequency 
has plateaus at 18 < n < 20, 24 < n < 26, 32 < ?i < 34, 
and then the deep minima. For large ?i ^ 31 -^- 32, the 
frequency becomes negative, lo <Q. The fact, that these 
negative values occur only for SIA for slightly less then 
the "magic" value n = 32, corroborates the above sug- 
gested concept. 

For the case of exchange anisotropy (Fig[33), the up- 
ward jumps on the curve uj{N) arc almost absent, and 
the smooth increase of uj occurs instead, i.e. the effect 
of "magic" numbers is absent. The downward jumps 
are clearly seen at the same values of n as those for 
SIA. These jumps sometimes are much sharper that those 
for SIA, but for ExA the amplitude of these jumps are 
smaller, than for SIA, and w > everywhere. 

The physical explanation of the above quite compli- 
cated behavior can be done based on the aforementioned 
picture. First, the upward jumps for n = rimag with 
subsequent almost constant u can be easily understood 
for SIA. In this case, at n = rimag the formation of fa- 
vorable coUinear structure has already been finished so 
that the plateau at larger n is due to the creation and 
growth of a "DW piece". It is clear that for ExA this 
scenario does not occur, and this effect is completely 
absent. The downward jump and general decrease of 
the soliton frequency is related to the transition to more 
symmetric configurations, where the excessive number, n 
of magnons is easily spread along the soliton boundary. 
Such transitions take place for both types of anisotropics, 
which explains the behavior similarities. 

Let us discuss now the possibility of the topological 
soliton existence in high anisotropy ferromagnets. In this 
specific case, it is useful to utilize a simplified obvious def- 
inition of the topological invariant. The 7r2 topological 
invariant l|12|) for the case of large anisotropy has simple 



FIG. 6: The same as in Fig. |^for larger values of ; 
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FIG. 7: Dependence uj{N) for a soliton with KeS = J, ob- 
tained from numerical simulations, (a)- Single ion anisotropy, 
(b) - exchange anisotropy. Dashed hues are for better com- 
parison and correspond to exchange anisotropy (a) and single 
ion one (b). 



geometric meaning. Namely, only 7^ 0, tt make a con- 
tribution into integral p2(l so that for its evaluation it is 
sufficient to consider the DW region. Formally, the inte- 
gral H12|) can be represented as a contour integral along 
the DW, 

This value defines a mapping of the DW line (which is 
of necessity a closed loop) onto the closed contour which 
is a domain of angle variation. This representation 
makes it obvious that the topological charge, Q notion is 
meaningful only in the case when the DW has a well de- 
fined noncoUinear structure throughout its length. The 
DW regions with coUinear spins play the role of a " weak 
link" , where </> can change abruptly by 27r almost with- 
out overcoming the potential barrier. Hence, even for 
nonsymmetric soliton textures, when soliton has a quite 
large "piece" of noncoUinear DW, with " quasitopologi- 



cal" spin inhomogencity, literally topological structures 
are absent. Although the difference between topological 
and nontopological solitons in the magnets with KeS ^ J 
is not that large and the question about realization of 
topological solitons in such structures is rather academic, 
this problem will be discussed in more details. 



32 
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2 0,3 OA 0,5 0,6 0,7 0,8 0,9 ^,0 K /J 

e 



FIG. 8: The energy of the soliton with N/S = 32 (magic) as 
a function of effective anisotropy constant. Inset shows the 
behavior of a derivative dE/dKcB near /('crit ~ 0.8. 

To answer the question about the presence or absence 
of a topological texture, we have carried out the nu- 
merical minimization both over a complete set of vari- 
ables and over 0's only with fixed in-plane spin direc- 
tions. Contrary to the above considered case of moder- 
ate anisotropy, the latter minimization (i.e. that over 0's 
only) never gives the instability of a nontrivial topological 
structure. Instead, the decrease of noncoUinear structure 
amplitude (either uniformly along the entire DW or on 
its individual parts) occurs as K^a increases. The be- 
havior is quite different for different values of n = N/ S 
so that these cases should be considered separately. For 
specific analysis, we have chosen a few n values, typical 
"magic" number n = 32 and two non-magic, n = 35 and 
n = 36. In spite of closeness of these numbers, the soliton 
spin texture behavior in these cases differs drastically as 
anisotropy increases. 

The energy of the topological soliton as a function 
of i^off / J for magic number n = 32 and two types of 
anisotropy is shown in FiglHl These curves have been 
obtained by numerical simulations on a lattice, for an 
anisotropy increase from its small value, when there is a 
well-defined topological soliton texture. It is seen that 
at large K^a the soliton energy tends to some finite lim- 
iting value, E'o = 32 JS"^, that is typical for a coUinear 
structure with the "magic number" n = 32. But the 
behavior of these functions near this limiting value is 
different for single- ion and exchange anisotropics, cor- 
responding to the different scenarios of annihilating of 
both non-coUinear spin structure and topological struc- 
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(a) ExA, 0.6 



(b) SIA, 0.6 




(c) ExA, 0.8 



(d) SIA, 0.8 




(e) ExA, 1.0 



(f) ExA, 1.08 



FIG. 9: Soliton structure for magic number of bound magnons 
N = 325 and different values of effective anisotropy constant 
Ke, the value of Ke/J, together with the type of anisotropy, 
are shown below pictures of spin distribution. 



ture in a soliton. It is interesting to trace the disap- 
pearance of a soliton topological structure at anisotropy 
increase, which is shown in Figs. El For small K^b the 
topological structure is well defined, an d this str ucture is 
similar for both SIA and ExA, see Figs.|9(^ 



) and 

the value of Kcs increases to some critical value 



9(bJ 



If 



the in-plane spin amplitude decreases, but a "vortex-like" 
configuration with approximate symmetry C4 is still vis- 
ible. For larger anisotropy K > A'crit the soliton struc- 
ture becomes purely coUincar. For SIA the critical value 
is ATcrit ~ 0.8, but for ExA this value is much larger, 
ATcrit ~ 1-1, compare Figs. 9(d) and 9(c) For ExA, 
the topological structure is still visible for s uch strong 
anisotropy as 2k = = J, sec Fig. 9(e) However, 



and of th e same structure as that for SIA, as is shown in 
Fig. |9(f)| In other words, for "magic" magnon numbers 
the topological structure of a soliton decays smoothly as 
the anisotropy constant grows, which is seen for both 



types of anisotropy, compare Figs. 9(d) and 9(f) 

For AToff > ATcrit for both kinds of anisotropy there is 
no more topological spin structure - all spins are directed 
either up or down and energy docs not depend on K^s 
any more (see FiglSJ. Such a picture resembles very much 
the phase transition of a second kind with the collinear 
soliton texture as more symmetric phase. This agrees 
with the behavior of the energy E{Kcff) near ATcrit, (see 
inset to FiglSJ 



(33) 



for high enough k the structure finally becomes collinear 



To understand better the picture of the transition from a 
topological (noncolhnear) soliton texture to a collinear 
texture, we discuss the analogy with the second or- 
der phase transition in more detail. Note, that any 
collinear spin structure with arbitrary positions of "up" 
and "down" spins has the same symmetry element, 
namely, rotation about the z - axis in a spin space. This 
symmetry is due to the symmetry of the Hamiltonian ^ . 
On the other hand, the spin textures, even noncolhnear, 
with n = 32 and C4 symmetry, are invariant with respect 
to a rotation by (7i'/2) • /c, fc g Z, simultaneously in spin 
space and coordinate space. Thus, "magic" collinear tex- 
tures with spatial symmetry C4 have higher symmetry, 
being invariant relative to the independent rotation of the 
coordinate space by (7r/2) • k and of the spin space by ar- 
bitrary angle, while noncolhnear "magic" ones are invari- 
ant only relatively to simultaneous rotation by (7r/2) • k in 
spin and coordinate spaces. This means that on a tran- 
sition from a "magic" collinear soliton to noncolhnear 
one, spontaneous symmetry breaking occurs and phase 
transition of the second kind appears naturally. 

For "non-magic" numbers 71 = 35 and n = 36, the 
symmetry of soliton structures and their behavior is fun- 
damentally different. First, note that the collinear struc- 
ture can be realized for even values of n only. For any 
non-even n, odd-integer or non-integer, some spins have 
to be inclined to the z— axis. Thus, these two cases have 
to be considered separately. 

The value n = 36 is rather far both from " magic" num- 
ber n = 32 and from the nearest "half-magic" n = 40, 
corresponding to the favorable configuration with a rect- 
angular collinear DW. Here, for small anisotropy K^s < 
0.6J we see the structure with C4 symmetry, similar to 
that for " mag ic" n ~ 32 for both types of anisotropy, see 
Fig. 9(a) and 9(b) However, with the increase of K^s, 
the evolution is different, as seen in Fig. ^| The dif- 
ference in behavior of E{Kcs) is the largest for exchange 
anisotropy, where the saturation, depicted in Fig. |H1 did 
not occur up to quite large k > 0.6J {KcS ^ 1.2J). How- 
ever, for single ion anisotropy there is also a difference 
from the above considered "magic" case n = 32. 

Those differences could be explained in terms of previ- 
ously discussed fact that DW pinning is stronger for SIA 
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FIG. 10: The energy of the soliton with N/S = 36 (non- 
magic, even) as a function of effective anisotropy constant 
Keff. For exchange anisotropy, fuU symbols presents data 
for topological solitons, open symbol non-topological solitons 
with Q — 0. The values of KcS « 1.13 for the transition from 
non-coUinear to coUinear soliton for single ion anisotropy are 
shown by vertical dotted lines. 
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(c) ExA, 0.1 



(d) SIA, 0.8 



compared to ExA. Using symmetry arguments, we note 
that the purely collinear texture for n = 36 has a DW of 



complex form, lower symmetry (Fig 11(h) I, and is much 
less energetically favorable, than that for n = 32. Hence, 
this state in its pure form occurs only for sufficiently large 
SIA with K > 1.18J rather then for K > 0.8J as for the 
"magic" case n = 32. We have found that the textures 
with sufficiently high symmetry (in any case they have a 
center of symmetry or a couple of C2 axes) are inherent 
to exchange anisotropy. Such textures for non-magic n's 
have DW's of strong non-coUinear structure with a great 
number of spins with 9 ~ 90°. For the ExA case, the 
DW does not "tear apart" up to quite large k > 0.5 J 
{KeS ^ J)- Here both topological and non-topological 
soliton textures exist. The energy of latter textures is 
shown in Fig. ^| by open symbols. For k < 0.4 both 
types of solitons have a rectangular shape and rather high 
symmetry. Their difference is that for non-topological 
soliton the inversion center is present even for KcS — 1.2 J 
while the symmetry of a topological soliton lowers al- 
ready for Kas — J, see Fig. ll(c)|ll(~ This difference 
can be understood as follows. For a topological soliton, 
the DW containing spins with 9 ~ 90° (since the spin 
angle </> is inhomogeneous) is less favorable than that for 
non-topological solitons. Thus, for topological solitons 
the spins with 9 ^ 0,tt concentrate near one of the soliton 
edges. However, the energy difference of these solitons at 
Kcff > I.IJ is quite small. Thus we may assert that 
in the wide range of ExA constants we have two soliton 
types with similar structure of spins along the easy axis 
and approximately equal energies. 

In the SIA case the DW pinning plays a much more 
important role so that the low symmetry, which is 




(e) ExA, 1.0 



(f) SIA, 1.0 




(g) ExA, 1.2 



(h) SIA, 1.2 



FIG. 11: Soliton structure for N=36 (even non-magic num- 
ber) and different values of effective anisotropy constant K^/ J 
(shown in captions) 



cha racteri stic of collinear texture of the type shown in 
can be traced already for sufficiently small 



Fig|TT(h) 



Kos > 0.7 J, see Figs |ll(d)|ll"(f]] Since the DW inho- 
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mogcneity, caused by a topology, is localized in a small 
part of a boundary, its role is not essential so that there is 
no big diflFerence between topological and non-topological 
solitons. At least this difference is much less pronounced 
than that for ExA where the non-coUinearity amplitude 
decreases smoothly as K increases. The new element of 
symmetry (a rotation in a spin space) appears in a tran- 
sition to the coUinear state at A'cs > ATcrit- That is why 
this transition resembles the phase transition of a second 
kind. This behavior of E{N), similar to that for "magic" 
magnon numbers, can be seen in Fig. llOl for SIA (but not 
for ExA!). 




FIG. 12: The energy of the soliton with N/S = 35 (non- 
magic, odd) as a function of effective anisotropy constant. 
Full symbols present data for topological solitons with inho- 
mogeneous DW, open symbols correspond to non-topological 
solitons with g = 0. 

As we have found out in the above example of n = 32 
and n = 36, the final stage of evolution of a soliton with 
even n is a coUinear structure. Since the transition of a 
soliton into a coUinear state is related to the symmetry 
increasing, it formally resembles a transition of the sec- 
ond kind. The coUinear structure cannot be realized, for 
obvious reasons, for any ri, which is not even. In this 
case, which has been considered for n = 35, the smooth 
variation of a soliton energy and structure occurs, see 
Fig. 1121 In particular, its energy is a monotonously in- 
creasing function of Kcft for both types of anisotropics. 
For high enough anisotropy, the difference between topo- 
logical and non-topological solitons diminish. 



VI. CONCLUSIONS 

In this paper we have studied the topological and non- 
topological solitons, stabilized by precessional spin dy- 
namics in a classical 2D ferromagnet with high easy-axis 
anisotropy on a square lattice. Our analysis has been 
performed both analytically in the continuous approxi- 
mation and numerically. The main conclusion is, that 



in the 2D case the solitons properties change drastically 
upon departure from the " singular point" - isotropic con- 
tinuous model with BP solitons. Similar to previous stud- 
ies, it turns out that the presence of even weak anisotropy 
makes solitons dynamic, i.e. those with nonzero preces- 
sion frequency for any number N of bound magnons. It 
is very interesting and unexpected, that the role of dis- 
creteness turns out to be of the same importance as that 
for anisotropy, even for k, K <ti 1, when /q S> a. The 
minimal consideration of discreteness (via higher degrees 
of gradients) yields the existence of some critical value 
(lower threshold) of both soliton energy and the number 
of bound magnons. Similar to the problem of cone state 
vortices, the instability is related to the joint action of 
discreteness and anisotropy "responsible" for Iq forma- 
tion (see Eq. (|14|l '). As a result, the nonanalytic, instead 
of linear, dependence of the critical soliton energy H31() 
on the anisotropy constant appears. 

On the other hand, for intermediate values of 
anisotropy, like < (0.25 0.3) J, and for the val- 
ues of N far enough from the critical value of the bound 
magnons A^cr, we have checked and confirmed a number 
of results regarding soliton structure from the continuous 
theory. In particular, the relation between the number 
of bound magnons and precession frequency of spins in- 
side the soliton is common to these two approaches. In 
agreement with Ref. ^^non-topological solitons have, as 
a rule, the lower energy than the topological ones (see 
Fig^, and they are stable for any N > Nnt- 

Some fundamentally new features in soliton behavior 
appears as the anisotropy grows. For "magic" numbers 
of bound magnons, N = 2/^5, I is an integer, the soliton 
texture corresponds to the most favorable structure with 
a coUinear DW of a square shape. In this case the soliton 
topological structure disappears (as anisotropy grows) 
continuously (to some extend analogously to the sec- 
ond order phase transition) giving purely coUinear state. 
Such behavior is inherent to both types of anisotropy, but 
for SIA the critical value Kait is lower and the transition 
is sharper. For non-magic numbers, when the coUinear 
structure is certainly less favorable, it still appears for 
even numbers N/S and the SIA case by a smooth tran- 
sition, but with essentially larger ATcrit- Note that in all 
cases this value of A'crit is much larger then that for ID 
spin chain, ATcrif = 0.5 J^^^. Finally, if n is not an even 
number, then for K > J the soliton always has a non- 
coUincar structure in the form of DW "piece" or even a 
single spin with Sz 7^ ±1- The energy of such a soliton 
is almost independent of inhomogencity of planar spin 
components in the DW. 

Let us note one more interesting property of the soli- 
tons in magnets with strong anisotropy, namely the ex- 
istence of the stable static solitons in this case. It is 
widely accepted, that such solitons are forbidden both in 
the model with W2 (Hobart-Derrick theorem) and in the 
generalized model with negative contribution of fourth 
powers of magnetization gradients {W2 + W4 in this arti- 
cle). However, the region of N values, where the soliton 
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frequency changes its sign, is clearly seen on the Fig. [7| 
for SIA. This means the existence of certain character- 
istic values of N, where uj = and the soliton is actu- 
ally static. Also, for given number of N in some region 
one can find the corresponding value of anisotropy con- 
stant (high enough) to fulfill the condition uj = and 
to realize this static soliton with non-coUinear state. In 
addition, all purely coUincar configurations, existing for 
sufficiently strong anisotropy of both types, are indeed 
static since the precession with the frequency lu around 
the easy axis does not define any actual magnetization 
dynamics. Hence, in the magnets with sufficiently strong 
anisotropy the stable static soliton textures can exist. 
But, contrary to Ref. they are due to discreteness 
effects, namely due to DW pinning. 

To finish the discussion of the problem considered in 
this work, let us consider the effects taking place at the 
transition from classical vectors Sn to quantum spins. 
Usually the simple condition of integrity of the total z- 
projection of spin A^, see ((T^ . is considered as a condition 
of semi-classical quantization of classical soliton solutions 
of Landau-Lifshitz equation for uniaxial fcrromagnetai. 
Moreover, the dependence E = ^^(A^) is used to be inter- 
preted as a semiclassical approximation to the quantum 
result. Note, that for some exactly integrable models, for 
example, XYZ— spin chain with spin S = 1/2, the ex- 
act quantum result coincides with that derived from the 
semi-classical approach, see Ref. 0). This simple picture 



was elaborated for solitons having the radial symmetry. 
Now we consider how this picture may be modified for 
solitons with lower spatial symmetry discussed in this 
article. Solitons with low enough symmetry, say "rectan- 
gular" solitons for "half-magic" numbers, or even lower 
symmetry, which occurs for " non-magic" numbers and es- 
pecially odd ones, apparently occur for high-anisotropy 
magnets. It is clear, that they can be oriented in a dif- 
ferent way in a lattice, that may be construed as fc-fold 
degeneracy of a corresponding state in the purely classi- 
cal case, with A: = 2 for a "rectangular" soliton or fc > 2 
for less symmetric states like those reported in the Fig. 
111! With the account for effects of coherent quantum 
tunneling transitions between these states, this degen- 
eration should be lifted. According to the semiclassical 
approach, valid for bound states of a large number of 
spin deviations in high anisotropy magnets, the transi- 
tion probability is low and can be calculated using in- 
stantons concept with Gaussian integration over all pos- 
sible instanton trajectories^. As a result one can expect 
the splitting of states degenerated in the purely classi- 
cal case, with creation of k-multiplet and lifting of the 
symmetry of the soliton. A detailed discussion of these 
effects is beyond the scope of the present work. 
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